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Abstract: We calculate the pion and p meson masses in quenched SU(iV) gauge theories 
for N = 2, 3, 4 and 6. Extrapolating these results to the chiral and large- iV limits, we find 
nip = (1.670 ±0.024)-y/cr for the p meson mass at a fixed lattice spacing a ~ 0.2093 a -1 / 2 ~ 
0.093 fm, where we use the (arbitrary) value (444 MeV) 2 for the string tension. We estimate 
a continuum limit large- N value, m p = (1.77 ± 0.05)^/a. We find 1/N 2 corrections to be 
small and we compare our results to predictions from AdS/QCD. 
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1. Introduction 

Numerical simulations of QCD in lattice regularization (Lattice QCD) are the standard 
way to derive low energy aspects of the phenomenology of strongly interacting elementary 
particles, directly from the QCD Lagrangian. Adding conformal invariance via supersym- 
metry or sending the number of colours N from three to infinity are popular ways of turning 
some aspects of non-perturbative QCD analytically tractable. Such extensions are in par- 
ticular also necessary to connect supergravity predictions to properties of induced four 
dimensional quantum field theories on the boundary of anti-de Sitter space |j|] (AdS/CFT 
correspondence or AdS/QCD @, || §). 

Pure SU(iV) gauge theories have been studied extensively on the lattice and, among 
other quantities, the glueball spectra, string tensions ||, MS A-parameters || and de- 
confinement transition temperatures |7| are well-established. These pure Yang- Mills results 
demonstrate that 1/N 2 corrections are already below the 10 % level at N = 3. Little is 
known as yet for full QCD with finite iV > 3, including np flavours of sea quarks, in which 
case corrections to the meson spectrum are expected to be only suppressed by factors 
n F /N. 
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Here, as a step in this direction, we study the meson spectrum in the quenched ap- 
proximation for N = 2, 3, 4 and 6. While these finite- iV results alone will not allow us to 
quantify the difference between N = 3 and the large- N limit of full QCD, the quenched 
and the full theory will share the same large- iV limit where sea quark loops are suppressed 
for any finite number of quark flavours ||] . 

Obtaining large- N results that are close to existing full N = 3 QCD results (or indeed 
to real world experiment) will help us to understand why the quenched approximation 
(and even the naive quark model) works remarkably well in many cases g. The large-iV 
expansion is also a powerful tool to reduce the number of low energy constants in chiral 
effective field theories 113] and we wish to understand in what cases three can be regarded 
as a large number and under what circumstances not. Moreover, QCD in the planar limit 
is an interesting quantum field theory in itself. Last but not least, only in this limit the 
AdS/QCD correspondence can work. 

First exploratory simulations of the pion masses in SU(17) and SU(19) were performed 
by Kiskis, Narayanan and Neuberger [11] some time ago, employing overlap fermions. 



Preliminary results for the tt and p meson masses using the Wilson gauge and fermion 
actions for N = 2,3,4 and 6 were presented by us at Lattice 07 [12]. The present article 



concludes this study. Recently, Del Debbio et al. |L3|, obtained the tt and p meson masses 
as well, using the same action at the same iV-values, at a somewhat coarser lattice spacing 
than ours. This coincidence gives us some control of the continuum limit. 

This article is organised as follows: in section ^ we introduce the methods employed in 
the simulation. In section |3] we describe our analysis methods and present the results. We 
discuss finite volume corrections, the ^-dependence of the additive quark mass renormal- 
ization and the p meson mass as a function of the tt mass at large N. In section |we relate 
our results to those of Ref. [|1| and perform the continuum limit extrapolation, before we 
conclude the article with a discussion in section ||. 

2. Simulation parameters and lattice methods 



We use the software package Chroma [14], to carry out our simulations which we have 
adapted to work for a general number of colours N. 
We use the Wilson plaquette gauge action, 

where U\j is the product of SU(iV) link matrices U X n, connecting the site x with the site 
x + ap, around an elementary plaquette □. a denotes the lattice spacing and p a unit 
vector in /^-direction. We use Wilson fermions with hopping parameter k, 

Sf = 5>jhk - «E " ijUx^z+afi + V4(l + l^Ul^A-afi] ■ (2.2) 

x X,fl 

k is related to the lattice quark mass m q by 

am --(-- — ) and k~ x = 8 + ©(/T 1 ) . (2.3) 

2 \K KrJ 
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To set the scale, we use the string tension calculations by Lucini et al. [15]. We choose 
the coupling [5 = 2N/g 2 = 27V 2 /A such that the infinite- volume string tension in lattice 
units ay/a is the same for each A. We use the value ay/a « 0.2093: for SU(3) this 
corresponds to f3 = 6.0175. A denotes the 't Hooft coupling in the lattice scheme. At 
our lattice spacing it converges towards the A — > oo value A 2.78. Adopting the value 
a = 1 GeV/fm ~ (444 MeV) 2 for the string tension, the lattice spacing is a ~ 0.093 fm in 
each case. The values of ay/a [[H| used in the fits are very accurate ||, so the mismatch 
of our lattice spacings in units of the string tension between different A is below the 1 % 
level. 

We have chosen our quark masses (or equivalently our values of k) in order to match 
the pion masses m n /y/a between the different A- values as closely as possible. These Re- 
values were estimated by carrying out short exploratory simulations of m n at different mass 
points and interpolating. 
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0.1457, 0.1480, 0.1500, 0.1510 


100 








24 3 x 32 


0.1510 


100 


3 


6.0175 


2.991 


12 3 x 32 


0.1547 


50 








16 3 x 32 


0.1500, 0.1520, 0.1537, 0.1547 


50 


4 


11.028 


2.902 


12 3 x 32 


0.15625 


50 








16 3 x 32 


0.1520, 0.1540, 0.1554, 0.15625 


50 


6 


25.452 


2.829 


12 3 x 32 


0.15715 


50 








16 3 x 32 


0.1525, 0.1550, 0.1563, 0.15715 


44 



Table 1: Simulation parameters. n con { denotes the number of independent gauge configurations. 



Our simulation parameters are given in table |. We carry out most of our calculations 
on lattices with volume 16 3 x 32 in lattice units, corresponding to a spatial extent of ~ 1.5 
fm in physical units. We carry out additional simulations on 12 3 x 32 lattices at our 
lightest quark masses, where finite size effects (FSE) are expected to be largest. In SU(2) 
we augment this by a 24 3 x 32 volume, since FSE are also expected to become larger with 
smaller A. Indeed, FSE are expected to be zero at infinite A, as long as the lattice is 



kept larger than a critical length l c [16]. Hence an extrapolation at a fixed finite volume 
to infinite N will yield the correct A = oo limit, in spite of any residual FSE we may 
encounter at finite A. 

We calculate correlators using three sources: 'point', 'narrow' and 'wide', and the 
corresponding three sinks. We employ Jacobi-Gauss smearing JO], 



(2.4) 



with smearing parameter 7 = 4, with 10 iterations for the 'narrow' sources and sinks and 
60 for the 'wide' ones. Prior to this we carry out 10 iterations of spatial APE smearing on 
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the parallel transporters appearing within eq. fl2.4|) above |l 



U'i = Proj 



SU(JV) 



4 Ux — a7,iUx 



(2.5) 



with smearing parameter a = 2.5, where U^{x) is projected back into SU(A) after each 
iteration. 

We invert the Dirac operator using the Conjugate Gradient (CG) algorithm with 
even/odd preconditioning. We find the number of CG iterations to be approximately 
independent of A at a given pion mass. Since each CG matrix-vector multiplication has 
a cost proportional to A and we have to solve for A different colour sources, the total 
cost of producing a propagator is proportional to iV 2 . Updating a gauge configuration 
(and also carrying out APE smearing) involves matrix-matrix multiplications and so has 
a cost proportional to ./V 3 ; however, for the relatively small values of A we use, the cost 
of calculating the propagators dominates, so the total computer time spent on generating 
and analyzing a fixed number of configurations is still roughly proportional to A 2 . 

We find that the correlators become less noisy as A increases, so we are able to extract 
masses with smaller errors from the same number of configurations. We find the errors for 
a given number of configurations to decrease approximately like 1/y/N (see tables [2] ||). 
Hence the number of configurations that would be required for equal precision decreases 
like 1/A and the total cost is only proportional to N. A heuristic argument supporting 
this observed reduction in noise is the increase in the number of degrees of freedom of the 
statistical system oc A 2 at fixed volume. 

However, in S it was observed that for glueballs both the signal and the error decrease 
as 1/A 2 , so the signal-to-noise ratio remains constant. Hence arguments based only on 
the degrees of freedom can be misleading and one must be more careful. Repeating the 
argument of ref. Q, one finds that a crucial difference is that when calculating the four-point 
correlator contributing to the noise the intermediate states are dominated by two-meson 
states, whereas in the case of glueballs the vacuum is a possible intermediate state. The 
upshot of this is that for mesons the error does decrease faster than the signal, by a factor of 
1/A, in agreement with the heuristic argument above. As noted above, the errors decrease 
by factors ~ 1/yN that are somewhat smaller than this best-case 1/A scenario. 



3. Results 

3.1 Extracting masses from correlation functions 

As described above, for each channel of interest we calculate cross-correlation matrices 
Cij{t) = (Oj(t)Oj(0)) ; where i, j correspond to the size ('point', 'narrow' or 'wide') of the 
sources and sinks. To extract masses, we first calculate eigenvectors ipf of the generalised 
eigenvalue problem M, [2(J : 

C-^o)C^o + a)Co=A«. (3.1) 

We then project C(t) onto the eigenvector ip a that corresponds to the largest eigenvalue 
A" and perform periodic cosh fits, varying the fit ranges, to extract the lowest mass. We 
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also vary to within the range - 2a; our results are not very sensitive with respect to the 
value of to- We estimate errors on the masses using the jackknife method. 
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Figure 1: Effective masses and fits for the tt ground states at the lightest mass- value for N = 2 
(k = 0.1510, m„ = 1.29(5)Vo : ) and N = 6 (k = 0.15715, m,, = 1.35(1)^- The lattice volumes 
were 16 3 x 32 in both cases and to = a for N — 2 while to = 2a for the more precise SU(6) data. 



As an example we show in figure |] the effective masses for our lightest pion in SU(2) 
and SU(6). We see that there are clear plateaus over a large range of t- values. These plots 
are typical of those for other N and other masses. 

As mentioned above (and as is also visible from figure |l]), we find that the correlators 
become less noisy, and hence the errors on the masses decrease, as N increases. We illustrate 
this in figure ^, where we compare point-point pseudoscalar correlators on individual gauge 
configurations in SU(3) and SU(6). The pion masses are almost identical, am^ = 0.276(4) 
and 0.283(2), respectively, but the scatter between individual configurations, a measure of 
the noise, is about twice as large in SU(3). This is in agreement with the naive degrees of 
freedom argument presented above. 

We have calculated correlators for every quark bilinear J PC channel. However, we 
find that only the J PC = h channel, corresponding to the pion, and the J PC = 1 
channel, corresponding to the p meson, have sufficiently strong signals to accurately extract 
masses at present statistics. We present the masses of the ground states in these channels 
in tables || [B|. 
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t/a 

Figure 2: Point-point pseudoscalar correlators on individual 16 3 x 32 gauge configurations at 
the lightest mass values in SU(3) (n = 0.1547, = 1.33(2) % /o : ) and in SU(6) (k = 0.15715, 
TOtt = 1.35(1)-^). The SU(6) results have been shifted vertically for the comparison. 
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0.279±J 
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Table 2: n and p masses in SU(2), in lattice units and in units of the infinite- volume string tension. 
The infinite volume results that are displayed for the pion in the last row are extrapolated. 

3.2 Finite size effects 

Finite volume corrections are expected to be largest at the smallest pion masses, and at 
the smallest values of N. Comparing our results for the lightest tt masses between different 
volumes (tables |2|-||), we see that the differences indeed appear to decrease with N. We 
also see that the corrections for the p meson are smaller than for the tt, as expected since 
the p is heavier. 

We can estimate the infinite-volume masses by fitting our finite volume results to the 



-6- 



K 


Volume 




/ / — 


am p 


/ / — 

mp/^a 


0.1500 


lb x 32 


0.552(3) 


o a a / o\ 

2.b4(2) 


n n c\r' ( c*\ 

0.60b( b) 


on ( \ 

2.89(3) 


0.1520 


16 3 x 32 


0.450(3) 


2.15(1) 


0.524( 5) 


2.50(2) 


u. ioo i 


ID X OZ 




1.66(1) 


0.448( 8) 


9 1 A ( A\ 
Z. 141 4 J 


n i ^47 

U. Id^ / 


1 9 3 x ^9 
1Z A OZ 


u.ouu^y j 


1.46(4) 


0.432(14) 


9 nfii'7 N i 




lb 3 x 39 

1U A. Oil 


?7ft(&\ 


1.32(2) 


0.398(10) 


1 90(5") 

l.t/UI J 1 




no 3 x 39 


D 974+ 4 

U.Z / *±_g 


1.3111 








Table 3: 


The same 


as table | for SU(3). 




K 


Volume 


am n 




am p 




n 1 r on 

0.1520 


lb x 32 


r\ r a a /o\ 

0.544(2) 


2.60(1) 


0.b00( 3) 


2.87(2) 


0.1540 


lb 3 x 32 


0.438(2) 


2.09(1) 


0.519( 4) 


2.48(2) 


n 

U. 1004i 


1 fi 3 v qo 


U.004i^O J 


1.69(1) 


0.459( 5) 


z. iy ^z j 


U. 100Z0 


1 9 3 x ^9 
1Z A OZ 




1.36(3) 


0.402(13) 


1 Q9l'fi"\ 

i.yz^o j 




lb 3 x ^9 

-LU A O^ 


u.zyoi o I 


1.41(1) 


0.424( 7) 


9 091^ 




no 3 x 39 


D 9Q4 +3 


1.40l 3 








Table 4: 


The same 


as table | for SU(4). 
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Table 5: The same as table | for SU(6). 



form, 

m^-(L) = m 7r (oo) [1 + C exp(— m 7r (oo)L)] , (3.2) 

with C positive |pl| , p2| . For SU(2) we have three values of L with which to carry out the fit, 
and we obtain am^ (00) = O.279I3 and C = 0.6±o;g. Note that the data can also be fitted 
to a constant. For the other values of N we only have two lattice sizes for a two-parameter 
fit, and furthermore the L = 12 lattices are rather small so higher order corrections to 
eq. ( [3.2D could be substantial. Hence for these cases we take the conservative approach 
1 of assuming that the coefficient C remains constant at its SU(2) value for all SU(iV). 
We have included the infinite-volume pion masses obtained from these fits in tables || |5l 
In no case do we obtain any statistically significant deviation between the infinite volume 
extrapolated result and those obtained on our L = 16a ~ 1.5 fm lattices. 

In principle we should also carry out similar corrections to the p meson masses. How- 
ever, the corrections are much smaller for the p, and will be negligible compared to our 

^ote that C -> as N -» 00. 
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statistical errors. So for the p we simply use the results from our largest volumes. Simi- 
larly, at larger masses it is not necessary to include finite-volume effects since they will be 
smaller than our statistical errors, for both the tt and the p. 

We note that our chiral extrapolations, described below, are not very sensitive to the 
details of our finite-volume corrections. This is because the chiral fits are mostly controlled 
by the data at higher mass values, which have significantly smaller errors, and are much 
less sensitive to potential finite-volume corrections. 

3.3 Determination of the critical hopping parameter 
3.3.1 k c at finite N 
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Figure 3: (am^) as a function of eq. (3.3) 



The Wilson action quark mass undergoes an additive renormalization (2clk c ) , see 



eq. d2.3|) . Thus we expect that the pion mass will be related to k (up to quenched chiral 
logs) by, 

(am w ) 2 



(3-3) 



By fitting to this equation we can extract k c for each N. We obtain good fits for each N, 
which we display in figure ||. Note also that the pion masses are horizontally aligned across 
the different iV-values, indicating that we have succeeded to approximately match the n- 
values to lines of constant physics, thus eliminating another possible source of systematic 
bias. The values of k c that we obtain are shown in table [6|. We note that while the lattice 
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0.05 0.1 0.15 0.2 0.25 

1/N 2 

Figure 4: n c as a function of 1/N 2 . For comparison we also show the infinite N extrapolation of 
the lattice 't Hooft coupling 2N 2 //3. 



't Hooft coupling 2N 2 / (3 at fixed string tension decreases with increasing N (see table [T]), 
the K c -values move away from the free field limit K c ,free = 0.125. 

We are dealing with the quenched theory and hence 
eq. (^^) should be modified by quenched chiral logarithms p3| 
giving 



[am-r, 



a- i 1 -!. 
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(3.4) 
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' 2 


0.1528(1)(13) 


3.0 
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0.1562(1)( 3) 
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0.1580(1)( 4) 


1.7 
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0.1588(1)( 2) 


1.7 


oo 


0.1596(2) 


0.3 



We also fitted our results to this parametrization but 
have found 5 to be consistent with zero, confirming that 
at m-x > 1.3y/a ~ 580 MeV we are not yet sensitive to 
quenched chiral logarithms. Comparing our best fits using 
eq. (|3.3j) to those where we allow 5 to vary, we estimate the 
uncertainty on k c due to the lack of control over the chi- 
ral logs to be significantly larger than our statistical errors 

(see table ^). This error should decrease with N since the quenched theory will become 
equivalent to the unquenched theory in the large-iV limit. Indeed, the N = 2 data allow 
for more curvature than the N > 2 data sets. 



Table 6: The critical hopping 
parameter k c , with the reduced 
X 2 of the fits. The first errors 
are statistical, the second sys- 
tematic. 



3.3.2 k c in the large- N limit 

We expect n c to have 0{1/N 2 ) correction to its large- iV value. Hence we fit the values of 
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table to the form, 



K C {N = oo) + ^ 



(3.5) 



After including the systematic uncertainties from the chiral extrapolation we obtain values 
ac c (oo) = 0.1596(2) and c = —0.028(3). Some of the systematics will be correlated and we 
obtain a rather small x 2 / n ~ 0.27. We display the 1/N 2 extrapolation of k c in figure ||. We 
also include the 1/N 2 extrapolation of the lattice 't Hooft coupling at fixed string tension. 
Note that this extrapolates to the value X(N = oo) = 2.780(4). 

3.4 The p meson mass 
3.4.1 m n at finite N 








2 3 4 



Figure 5: Fits of m p /^fa as functions of rn^/a, eq. (|3.6[), at different N 



In chiral perturbation theory, as well as in 
the heavy quark limit, m p depends linearly on the 
quark mass m q . Within our range of pion masses 
"T-vr/v^ 7 ~ 1-3... 2. 6 we find m 2 , to linearly de- 
pend on ft -1 and therefore to be proportional to 
the quark mass. Hence, we can fit our p masses to 
the parametrization, 
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1.64(3) 
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0.6 



Table 7: The fit parameters of eq. (3.6) 
for each N, with the reduced x 2 -values. 
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where m p (0) denotes the p meson mass in the chiral limit. Quenched chiral logarithms will 
modify the relationship between m q and m 2 , however eq. ( |3.G| ) remains valid |23( |. 

We obtain good fits in each case, which we show in figure |H[ We display corresponding 
parameter values m p (fi) / yfo and B in table [/]. 

3.4.2 m p in the large-iV limit 
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Figure 6: m p /y/a in the chiral limit, as a function of 1/N 2 . 

In the quenched approximation the parameters m p (0) / ' ^fa and B are each expected 
to have 0(1/N 2 ) corrections to their large-iV values, in analogy to eq. (|3.5|) . We see from 
figure H that the results and fits for different ./V are all very similar, with all points nearly 
lying on one line, so these corrections should be small. We obtain, 



m p (0)/y/a = 1.670(24) 



0.22(23) 



N 2 



and 



B = 0.182(5) 



0.01(5) 

N 2 



(3.7) 



(3.8) 



The fits are both good, with reduced x 2 -values of 1.8 and 1.3, respectively. We plot the fit 
for m-p(O) / ' \fa in figure |6|. Note that the iV-dependence is extremely small both for m p /yfa 
and for B. In particular, the difference between m p (0) for SU(3) and SU(oo) amounts to 
only (1.4 ± 1.6) %. 



Equating ^Jo 444 MeV, the large- N chiral p mass, eq. (3.7), corresponds to 



m p {m w = 0, N = oo) » 741(11) MeV 



(3.9) 
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The above value a ~ lGeV/fm is motivated by Regge trajectories, potential model fits 



to quarkonium spectra and unquenched N = 3 lattice data [24|. However, there is a 
systematic uncertainty associated with it, in particular also since the QCD string at finite 
N and np > will decay, once a critical distance is reached [pq] . Moreover, no N = oo 
experiment exists. Nonetheless, we find it interesting that this mass-prediction for the 
stable N = oo p meson comes close to the mass m p ~ 775 MeV of the experimental p 
resonance, which in fact has a significant decay width, V ~ 150 MeV. Note that at physical 
7r mass we obtain the large- N limit m p ~ 749(11) MeV. The (unknown) scale uncertainty 
from the string tension is not included in the above error estimates and neither are finite 
lattice spacing effects. The continuum limit extrapolation of section |] below increases our 
large-iV m p prediction by another ~ 45 MeV. 

3.5 Excited states 

In principle, our correlation matrices enable us to extract the masses of excited states as 
well as of the ground states. However, in practice we obtain poor mass plateaus for the 
excited states, and our statistical errors are rather large. Having said that, our results are 
consistent with the first excited pseudoscalars and vectors having no iV-dependence at the 
10 % level. Also, we find that our correlators at a given pion mass are similar for all N, 
not just at large t where they are dominated by the ground state, but all the way back to 
t = where many excited states contribute. An example of this can be seen in figure |2|. 
This is consistent with the entire spectrum of excited states to only weakly depend on N. 
The correlators also depend on the overlaps of our operators with the excited states, and 
these overlaps in turn depend on wavefunctions of the excited states. So this also suggests 
that the meson sizes and internal properties do not change strongly with N. 

4. The continuum limit 

Presumably there are lattice spacing corrections to our results. We can get some idea of 



their size from the study of Del Debbio et al. [13], on somewhat coarser lattices. 

Turning first to the large- N limit of their results, and using the facts that their lattice 
spacing is given by a as 1/(5 T c ), where T c denotes the de-confinement temperature, and 
that at A^ = oo one has T c /J~a w 0.5970 [15|, they obtain, 

2 

777 777 

-£ = 1.609(9) +0.1750(3)— ^ ; 0^ = 0.3350. (4.1) 
Jo u 



We can compare this to our large-A" result, eqs. (|3.6j) -(|3^8| 



—L = 1.670(24) + 0.182(5)— ^ ; o^ff = 0.2093 . (4.2) 
'a o 



We see that, despite the 60 % difference in lattice spacings, the results for both the chiral 
limit and the slope of the p mass are very similar, differing by less than 5 %. Also the ratio 
of our and their results remains remarkably independent of the ir mass, when expressed in 
units of a. 
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The leading lattice artefacts are of order a and a simple linear extrapolation yields, 

,2 

; (a - 0) . (4.3) 



mz 



1.77(5) +0.193(14)- 
a a 



Obviously, with just two values of the lattice spacing, we have little control over the system- 
atics of this continuum limit extrapolation. Therefore we also attempt a purely quadratic 
extrapolation in a and find the central values 1.70 and 0.186 for the two above parameters, 
respectively. Based on this and results of previous large scale N = 3 quenched spectroscopy 
studies with Wilson action Q , that included similar lattice spacings (as well as finer ones) , 
we conclude that the statistical errors stated above are sufficiently large to accommodate 
the possible contributions from subleading terms. 
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Figure 7: Results for m p /yfa and rn^/a from this work (full symbols, a sa 0.209/ y/a) and from 
ref. Jl3| (open symbols, a = 0.2/T c as 0.335/ \/o). The lines correspond to the respective N — > oo 



limits and the grey error band denotes our large- TV continuum limit estimate eq. (4.3) 



We now turn to the finite- N corrections. Here the two sets of results appear to be rather 
different. Our results show a very weak iV-dependence, with for example the difference 
between m p (0) for SU(3) and SU(oo) being only (1.4 ± 1.6) %. The corresponding value 
of ref. |L3| is (12.8 ± 0.6) %, much larger. Most of this difference can be attributed to the 



different ways of setting the scale; the ratio T c /y/a depends on N |15]. To remove this 
effect we have rescaled the results of ref. jnj, using the appropriate values of T c /^fa for 
each N [15|. We compile all these results in figure [7|. Again we remark that the results 
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of ref. were obtained at different values of ay/a for the different A. Also note that 
the range of pion masses covered in this case varies with A. Nevertheless, we see that 
when plotted in this way the two studies agree very closely, with the exception of SU(2). 
This suggests that most of the difference between the two data sets is due to the different 
methods of setting the scale, rather than from lattice corrections. The finite lattice spacing 
effects might be larger for SU(2), accounting for the deviation there. 

The fact that the iV-dependence at our finer lattices is already consistent with zero 
leaves little room for 1/A 2 corrections to m p /y/a being larger than 3 % at A = 3, at least 
for pion masses m n / ' y/a < 2.6. The figure also contains the large-A extrapolations of the 
data sets at the two lattice spacings and our continuum limit estimate eq. ( fO|) (grey error 
band). 

All this suggests that our finite-a results deviate by no more than 5-10 % from the 
continuum limit. This is true both for the large-A limit, eq. ( |4.2| ), and for the finite- A cor- 
rections. Also, these corrections are smaller when the masses are expressed in terms of the 
string tension than when they are expressed in terms of the de-confinement temperature. 
Note that this is opposite to the situation with respect to the scalar glueball which scales 
better between different A when normalized with respect to T c , rather than by y/a ||. 



5. Summary and Discussion 

We have studied the p and it meson masses in quenched SU(A) QCD with A = 2, 3, 4, 6 
at one value of the lattice spacing ay/a = 0.2093, in units of the string tension, at four 
different values of the quark mass. Our main result can be summarized by combining 
eqs. (pi), <K% and (O), giving 



"VKr) _ i °' 22 ( 23 ) , Lso/fl °- 01 ( 5 ) 



1.670(24) - + 0.182(5) - -f . (5.1) 



m 



2 



Combining our data with that of Del Debbio et al. allows us to extrapolate the large- A 



result to the continuum limit, see eq. (^4.3|) . 

We can compare our results to predictions from AdS/QCD correspondence. As an 
example, we consider the case of the Constable-Myers deformation analyzed in ref. 0|. 



In their original units of m p (0) this reads [27] 



mJiriTr) ( rn-n , 
1 " - 1 + 0.307 £-r . (5.2) 



m p (0) \m p (0) 
At our value of a we obtain 0.304(3) for this slope in units of m p (0), while the values 



obtained in ref. [13] can be translated into a coefficient 0.2816(2). A linear continuum limit 
extrapolation of these two numbers yields 0.341(4) and a quadratic extrapolation results 
in 0.318(3). The expected dominant behaviour is linear and hence we quote 0.341 ± 0.023 
as our large- A continuum limit result, where the error is dominated by the systematics 
of the continuum limit extrapolation. The above AdS/QCD prediction eq. ( [5, 2D agrees 
reasonably well with this result. 
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In the large- N limit the quenched theory is equivalent to the full, unquenched, theory. 
We have found that the 0(1/N 2 ) quenched corrections to the large- N result are very 
small. What about the 0(np/N) corrections of the unquenched theory? We can get an 
idea of their size without needing to carry out full unquenched simulations, by comparing 
our results with experiment. Converting the a — * and N — > oo extrapolated result of 
eq. ( [4 .3[ ) into physical units gives, 

9 

TTh 

m p (m„) = (786 ± 22) MeV + (435 ± 34)—^- . (5.3) 

GeV 

At physical tt mass this yields an infinite N value, m p = (794 ± 23) MeV, in perfect 
agreement with experiment: m p ~ 775 MeV. This might suggest that not only the quenched 
l/iV 2 -corrections but also the unquenched rcp/iV-corrections are small at N = 3. 

Of course the string tension value sfa = 444 MeV is arbitrary and we did not associate 
any error to this choice. Still this is very encouraging and might explain the success of the 
quenched approximation in light hadron spectroscopy calculations of the 1980s and 90s ||. 
It is certainly worthwhile to extend the present study to other mesonic channels, to confirm 
this picture. 

As discussed above, the fact that our results agree reasonably well with those of ref. [13] 



although our lattice spacings differ by as much as 60 % suggests that the systematic error 
due to the finite lattice spacing is already small at a ~ 0.093 fm for the Wilson action. 
Our continuum limit extrapolation confirms it to be smaller than 10 %. It would be nice 
to further constrain the continuum limit by repeating our work with the same analysis 
methods at another, smaller, lattice spacing. 

The errors associated with the neglect of quenched chiral logs in our pion mass fits could 
be reduced by going to lighter tt masses; the obstacle of exceptional configurations should 



be reduced at high N since the Dirac matrix eigenvalue distribution becomes narrower [12]. 
Our calculations have been carried out on more than one volume, so we already have some 
control over the small finite volume corrections. These should in any case disappear in the 
large- iV limit. 

We have not been able to obtain accurate masses for excited states in this work, 
although our results do suggest that they do not have a strong iV-dependence. We intend 
to use improved meson operators to calculate the masses of the low-lying excited states 
more accurately in the future. This should yield several independent mass ratios in the 
large- iV limit which can be used to constrain AdS/QCD models. This would also shed 
more light on the question if the use of the quenched approximation at np/N = 3/3 is 
justifiable from a large- N viewpoint and relate low energy constants in chiral Lagrangians. 
We are also planning to study the scalar sector and flavour singlet diagrams. 
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